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Abstract 

Nonperturbative dynamics of quantum fields out of equilibrium 
is often described by the time evolution of a hierarchy of correlation 
functions, using approximation methods such as Hartree, large Nf, 
and nPI-effective action techniques. These truncation schemes can 
be implemented equally well in a classical statistical system, where 
results can be tested by comparison with the complete nonlinear evo- 
lution obtained by numerical methods. For a 1 + 1 dimensional scalar 
field we find that the early-time behaviour is reproduced qualitatively 
by the Hartree dynamics. The inclusion of direct scattering improves 
this to the quantitative level. We show that the emergence of nonther- 
mal temperature profiles at intermediate times can be understood in 
terms of the fixed points of the evolution equations in the Hartree ap- 
proximation. The form of the profile depends explicitly on the initial 
ensemble. While the truncated evolution equations do not seem to 
be able to get away from the fixed point, the full nonlinear evolution 
shows thermalization with a (surprisingly) slow relaxation. 
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1 Introduction 



An understanding of the dynamical evolution of nonequilibrium quantum fields 
is needed in diverse physical situations as nonrelativistic condensed matter, rel- 
ativistic heavy-ion colliders, and the early universe. For several reasons a theo- 
retical description is expected to be difficult. First, effective irreversibility has 
to arise from time-reversal invariant equations. Second, for (asymptotically) late 
times one expects a nonequilibrium system to thermalize, which implies an ef- 
fective independence of the initial state, and prescribes a definite value for all 
correlation functions. Furthermore, during the nonequilibrium evolution, the an- 
swer to the question: what dominates the full dynamics at a certain stage, might 
itself be time-dependent. When approximation methods are used, the chosen 
method will often need to be modified or replaced as time goes on, to incorporate 
this shift in importance. 

To illuminate this, consider as a prototype for a nonequilibrium system the 
universe at the end of inflation (see e.g. 0,3)- In this case several distinct 
regimes are easily identified. While the early stage is dominated by the oscillating 
inflaton field, in the intermediate regime the created quanta scatter both with the 
inflaton and with each other, leading to a partial energy redistribution between 
the modes. Finally, in the last stage subsequent interactions are expected to 
bring the universe to thermal equilibrium. It is a great theoretical challenge to 
describe these various stages in a unified way. 

Although in principle the time evolution of expectation values is determined 
completely, after the initial density matrix is given, by the microscopic Heisenberg 
equations of motion, this is in practice only of minor help, due to the absence 
of exact solutions or solution methods. Consider for example a simple scalar 
quantum field theory with a quartic interaction. The Heisenberg equation of 
motion determines the time evolution of the mean field, 

{d 2 t -V 2 X + m 2 ) (<f>(x, t)) = -A(0 3 (x, t))/2, (1) 

where the brackets indicate the expectation value with respect to the initial 
density matrix. This equation requires knowledge of the three-point function 
{4> 3 (x,t)). The equation for the three-point function itself involves the five-point 
function (in general n + 2-point functions are needed to solve the exact evolution 
for n-point functions), which leads to a full hierarchy of coupled equations. In 
most solution to this hierarchy is not available, and it becomes of interest 

to find approximation schemes that capture the essential part of the full dynamics 
as correctly as possible. 

A widely-used approach is to truncate the infinite hierarchy of correlation 
functions. One of the simplest truncations is the Hartree approximation in which 
at most two-point functions appear, and the three-point function is replaced by 
(4> 3 (x, t)) = 3((j)(x, t))(<j) 2 (x, £)). A systematic way to implement this is by using a 
large Nf expansion, where Nf denotes the number of e.g. scalar or fermion fields 
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H (for applications, see e.g. ||, ^ || 0). The main drawback of these Gaus- 
sian truncations is that in homogeneous or translationally invariant ensembles 
scattering is absent, which limits the range of validity. This feature has been an 
important stimulation to improve upon the homogeneous Hartree and large Nf 
approximations. One possibility is to allow for inhomogeneous mean fields. In 
that case scattering via the space-dependent mean field ((f>(x,t)) is present in the 
effective equations (see e.g. || for analytical investigations and |§ for a numerical 
study). Another natural extension is to go beyond the Gaussian approximation 
and include higher nontrivial correlation functions |L0|, |Tl|, [l^], guided e.g. by the 
large Nf expansion to next-to-leading order. This typically results in effective 
equations that are nonlocal in time. From a practical (numerical) point of view, 
it is desirable to use effective equations that are local in time, as in the Hartree 
approximation. This has motivated the use of the lPI-effective action for equal- 
time correlation functions 13, I3|. Truncated at quartic order this includes 



scattering and all 1/Nf corrections fl6| . A discussion of the structure of dynam- 
ical equations for equal-time correlation functions in the large Nf expansion up 
to and including 1/Nj terms can be found in |fLT|| . Finally, very promising results 
have been obtained recently [|TI|, using a truncation of the 2PI-effective action 
at next-to-leading order. The resulting equations are nonlocal in time, i.e. they 
require the integration of memory kernels. For a 1 + 1 dimensional field theory 
this appears, however, not to be a numerical obstacle. 

In all cases, it would be desirable to perform a direct test of these methods and 
their range of validity. Unfortunately a comparison with the full nonperturbative 
evolution in the quantum field theory cannot be made, due to the absence of 
exact methods. In the case of quantum mechanics on the other hand, such 
tests can be performed, and the evolution from a Hartree factorization, a 1/Nf 
expansion at leading and next-to-leading order, and other extensions beyond the 
Gaussian aproximation have been compared with the exact evolution obtained by 



numerically solving the Schrodinger equation [|10| |TT[| . It is, however, not obvious 
how the lessons learned from quantum mechanics with one degree of freedom 
can be translated to field theory with (in principle infinitely) many degrees of 
freedom. Both scattering and the possibility of taking the thermodynamic limit 
are absent in the quantum mechanical case. 

The situation in classical field theory is completely different. Here the full 
evolution can be simulated using Monte Carlo methods and numerical integration: 
initial conditions are generated by sampling according to the initial probability 
distribution, and the subsequent time evolution follows from solving the classical 
equations of motion, which can be done numerically. Expectation values are then 
constructed by summing over many independent realizations. When the number 
of initial conditions is taken larger and larger, the initial probability distribution 
is approximated better and better, and the resulting time evolution will become 
exact, in principle. As we will see below, it is possible to implement many of the 
approximation methods discussed above for the quantum field theoretical case 
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also in a classical field theory. The reason is that the methods do not directly 
touch upon the quantum nature, but instead state how to truncate a hierarchy 
of correlation functions, which is present in a classical statistical system as well. 
Therefore, we focus in this paper on nonequilibrium evolution in classical field 
theory, formulated on a spatial lattice to regularize the theory. Note that the role 
of the thermodynamic limit may be investigated keeping the lattice spacing fixed. 
The general strategy is to compare the truncated dynamics with the numerical 
results obtained from sampling initial conditions from a given initial probability 
distribution. Our hope is that the insights obtained below will survive when 
going to the quantum field theory. 

In the remainder of the Introduction we give the outline of the paper. As a 
simple toy model we consider a classical scalar field theory in 1 + 1 dimensions, 
with the (continuum) action 



In the next section we show how the (homogeneous) Hartree approximation can 
be implemented in the classical theory. In Sec. [3] we discuss evolution equations 
obtained from a functional differential equation for the equal-time 1PI effective 
action. A truncation of the time-dependent effective action at quadratic order 
gives evolution equations that are equivalent to those obtained in the Hartree 
approximation. Scattering is incorporated by a truncation of the effective action 
that includes momentum-dependent four-point functions. In Sec. (| the choice 
of initial probability distribution and some numerical and lattice aspects are 
discussed. We have divided the comparison between the truncated evolution 
equations and the lattice results in two parts. In Sec. |5] we discuss the early 
and intermediate times, where the system is still relatively far from thermal 
equilibrium. The late-time regime, where thermal equilibrium is approached, is 
described in Sec. ||. Our findings are summarized in Sec. [7|. 

2 Classical Hartree approximation 

The classical problem is fully specified by the equation of motion 



supplemented with initial conditions for (j)(x, t) and 7t(x, t) = dt<p(x, t). These ini- 
tial conditions are determined by the initial probability distribution p[ir(x), 4>(x)}, 
where (p(x) = (fi(x,0),7i(x) = 7r(a;,0). The choice of initial distribution is not 
needed at this stage, we only assume that the ensemble is translationally invari- 
ant in space and respects the discrete symmetry <fi — > — 0, 7r — > —n. The average 
with respect to the initial distribution will be denoted with brackets (•). Our aim 





d 2 t (f){x,t)= d 2 x -m 2 <j>(x,t) - \<p 3 (x,t)/2 



(3) 
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is to find evolution equations for (equal-time) correlation functions of the field 
4>(x,t) and the canonical momentum ir(x,t). 

As discussed in the Introduction, a widely-used approach in quantum field 
theory is the Hartree approximation. This can be implemented in the classical 
theory as well, and a brute but simple way to do this is as follows. If we assume a 
Hartree-type factorization for the interaction term in the classical equation, i.e., 
\(j) 3 -»■ 3A0(0 2 ), we find 



d 2 <f)(x, t) 



dl-m 2 -\\^{x,t)) 



(x,t). (4) 



In the case of translationally invariant ensembles, (<p 2 (x,t)) is independent of x 
and the equation can be written in momentum space as 

%<t>(q,t) = -Qfy(q,t), (5) 
with the effective frequency squared 

-? = ^ 2 + ^(0 2 ), (6) 

where uj 2 = q 2 + m 2 , and (<ft 2 ) = (<ft 2 (x,t)). In the case of Nf scalar fields with 
a complete O(Nf) symmetry, the second term on the right-hand-side of Eq. (|6|) 
is multiplied by (Nf + 2)/(3Nj). In this paper we restrict ourselves to the case 
Nf = 1. Note, however, that there is no problem in principle to extend the 
analysis below to finite Nf > 1. The unequal-time two-point function 

S(x - y; t, = t')) = [ ^ e^^S(q- 1, t'), (7) 

obeys in this approximation the usual mean-field equation of motion 

d 2 +LJ 2 ]S( q] t,t') = 0, (8) 

where the effective frequency is determined from the two-point function at equal 
time: 

u 2 q =u 2 q + \\j^S( V] t,t). (9) 

Note that modes with a given momentum q only interact with the homogeneous 
mean-field background, so that direct scattering between different momentum 
modes is absent. 

The dynamics can be written equivalentlyQ in terms of equal-time expectation 
values, at the expense of introducing more than one two-point function. The 



the sense that the evolution equations presented below are obtained from the same 
starting point (||). 
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following four combinations are a priori independent: 
G<M>{x -y,t)= ((j)(x, t)<j)(y, t)), 

G^{x -y,t) = (n(x,t)n(y,t)), (10) 
G^{x -y,t) = -(rc(x, t)<f>(y, t) + <j>(x, t)n(y, t)), 

and the parity-odd combination 

G°f(x -y,t) = (n(x, t)<f>(y, t) - <f>(x, t)n(y, t)). (11) 

The dynamical equations are conveniently written in momentum space, according 
to 

G H {q,t) = J dxe~ iqx G^(x,t), (12) 
etc., and they couple only the combinations (110): 



d t G^(q,t) = 2G n ^(q,t), 

d t G W(/) (q, t) = -CjgG^q, t) + G n7T (q, t), (13) 
d t G^(q,t) = -2u^G^{q,t). 

The fourth combination G°^(q) is exactly conserved under the Hartree equations 
and does not enter in the dynamics. It is zero in the case that the ensemble is 
invariant under space reflection. 

In the Hartree approximation a nontrivial combination of two-point functions, 
termed a 2 in | 14| , is conserved, 

a-\q) = G^(q, t)G^(q, t) - G^q, t), (14) 

for each q. This can be understood from the absence of scattering or mode mixing 
at this order and the resulting symmetry (see Appendix |A|). Finally, the Hartree 
equations conserve the expectation value of the energy, 



f dq [l I/2 3A f dp , \ , 

^Hartree = L J — -G nw {q, t) + - [ UJ q + — J — G^(p, t) I G H {q, t) 



, (15) 



which can be obtained using a Gaussian factorization of the 4 -term in the mi- 
croscopic expression for the energy, or from the effective Lagrangian given in 
Appendix [A]. 



3 Equal-time effective action 

In order to improve upon the Hartree approximation, it is necessary to include 
direct scattering contributions in the evolution equations. We aim here at a 
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formulation that is local in time. This is also desirable from a numerical point 
of view.0 One way to achieve this is by employing a formalism based on the 
equal-time effective action T[<f>, 7r;t], the generating functional of 1PI equal-time 
correlation functions [13|. The effective action obeys the following (functional) 

(16) 



evolution equation [13, [16 1 



d t T[(f),ir;t] = -£ cl r[0,7r;t], 



with 

£ C 1 



dx 



ir(x) 
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X 



+ 4>{x) 



dl-m 2 - L -X 
x 2 



(j) 2 {x) + SG^x, x) 



5 



5n(x) 



dx±dx 2 dx 3 G<wi xifG^ (x, x 2 )G^ 3 (x, x 3 ) 
5 3 T 5 



Sip 1 (x 1 )5ip 2 (x 2 ) Sip 3 (x 3 ) 5ir(x) 



;i7) 



Here if) = (0, tt) and G^(x, x') denotes the full (matrix) propagator in arbitrary 
field background, obtained from V as 



6 2 T 



6ip(x)ip'(y)' 



(18) 



The full propagator evaluated at zero background is written without the bar. The 
effective action depends on 'effective' fields and tt, which are defined via the 
Legendre transformation, relating T and InZ in the usual way |13|]. Though we 
use the same notation, these fields should not be confused with the microscopic 
fields that appear in the original action (fj). 

In order to solve the exact evolution equation (|16|), some approximation has 
to be made. This brings us back to the issue of truncations, as discussed in the 
previous sections. We use a truncation or ansatz that includes all 1PI n-point 
functions, with n < 4, and respects the symmetry if) — > —if) as well as spatial 
translation and reflection. The ansatz reads 



r[0,7r;t] 



^(g)0*(g)0(g) + ±B(q)n*(q)n(q) + C(g)vr*(g)0(g) 



1 



2ir5(q 1 + q 2 + q3 + qi) u(q x , q 2 , g3)0(<?i)0(g2)0(<?3)0(<?4) 

51>92,?3,?4 

+v(qi, q 2 , q3)n(qi)(f)(q 2 )(f)(q3)<i)(q4) + w(q u q 2 , gaM^iM^)^)^) 
+V(qi, <12, q?>)K{qi)K(q 2 )Tx{q 3 )<f){qi) + z(q u q 2 , gsMgiM^MgaM^) 

We use the shorthand 

dq 
2^' 



2 See however uM for a successful implementation of time-nonlocal evolution equations. 
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and suppress the time dependence of the two-point couplings A,B,C, the four- 
point couplings u, v, w, y, z, and the correlation functions in this section. 
The matrix relation (pl|) at vanishing background can now be given explicitly, 



and 



G H {q) = B(q)/a 2 , 

G^(q) = -C(q)/a 2 , (19) 
G^(q) = A(q)/a 



2 



a 2 {q)=A{q)B{q)-C 2 {q). (20) 



with the determinant 



This definition of a(q) is equivalent to that in Eq. ([14]). For future convenience 
we introduce 

which will be used to convert two-point functions: 

G^(q) = -c{q)G^{q). (22) 



The time dependence of the effective action determined by ( |16|) translates 
into evolution equations for the couplings. Exact flow equations for the two- 
point functions follow from taking the second derivatives of Eq. ( |16|) with respect 
to and 7r at = 7r = 0:Q 

d t A(q) = 2Co 2 C(q) 

d t B{q) = -2C{q) - 2 1 {q)B{q) (23) 
d t C{q) = -A(q) + Co 2 B(q) - >y(q)C(q), 

with the frequency squared 
3A 



^ + yI p gm (24) 

3A f 

— 5" / 2n5 ( q _ ?i _ ?2 - G^{q l )G^{q 2 )G H {q z ) Au{q x , q 2 , q 3 ) 

O Jq\,q2m 

-?>c(qi)v{q u q 2 , q 3 ) + 2c(q 1 )c(q 2 )w(q 1 , q 2 , q 3 ) - c(q 1 )c(q 2 )c(q 3 )y(q 1 , q 2 , g 3 )] , 

where we recognize the free part, the Hartree term, and a contribution with 
the topology of the setting-sun diagram, containing three full propagators and a 



3 The equations presented below are slightly simpler than the ones that can be obtained from 
setting Nf = 1 in the equations for the O(Nf) model |16fl. 
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COq = ( 



Figure 1: Graphical representation of the effective frequency squared, Eq. 
The thick (thin) lines denote the full (free) equal-time two-point function G^. 
The blob in the setting-sun diagram is a full equal-time vertex function. 

(complicated) dynamical vertex function (see Fig. p]). The other factor appearing 
in ( P3| ) has a setting-sun structure as well: 

3A f 

= T 2 ^^(q- q t - q 2 - q 3 )G^(q 1 )G H (q 2 )G cf>4> (q 3 ) v(q,-q 1 ,-q 2 ) 

O Jq 1 ,q 2 ,q3 

-2c(gi)w(-gi, g, -g 2 ) + c(g 2 )c(g 3 ){?/(g 3 , g 2 , -g) + 2y(-g, q h g 3 )} 
-4c(gi)c(g 2 )c(g 3 )z(-g 1 , g, -g 2 )] . (25) 

These equations have to be completed with the evolution equations for the four- 
point couplings, and those are listed in Appendix [B|. Note that a(g) is no longer 
a conserved quantity as it is in the Hartree approximation, but obeys d t a(q) = 
-j(q)a(q). 

Finally, the evolution equations conserve exactly the energy 
E r = -^Hartree - / 2n5(q 1 + q 2 + q 3 + q±) G^(gi)G00(g 2 ) 

xG , ^(g 3 )G^(g 4 ) ]u{qi, g 2 , q 3 ) - c(gi)f (gi, g 2 , g 3 ) + c(gi)c(g 2 )w(gi, g 2 , g 3 ) 

-c(gi)c(g 2 )c(g 3 )?/(gi, g 2 , g 3 ) + c(g 1 )c(g 2 )c(g 3 )c(g 4 )2;(g 1 , g 2 , g 3 )] . (26) 

It is illuminating to make a connection with the Hartree equations derived in 
the preceding section. A truncation of the effective action at quadratic order gives 
dynamical equations involving only A, B, and C. It is straightforward to check 
that these give precisely the Hartree equations (0) for the equal-time two-point 
functions. 

The equal-time effective action permits an easy inclusion of quantum effects 



]j| . The Hartree approximation does not distinguish between classical and quan- 
tum field theories. In the quartic truncation the quantum effects add simple terms 
to the evolution equations of the quartic couplings |14[]. A direct verification of 
the truncated evolution for quantum fields is obviously much harder. 



4 Initial ensemble and lattice discretization 

For an investigation of the time evolution of correlation functions we need to 
specify the initial ensemble or probability distribution. In this paper we choose 
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to start from a Gaussian, translationally invariant ensemble. The reasons to 
take an initially Gaussian ensemble are the following: first of all, the Hartree 
equations truncate the dynamics to Gaussian dynamics for all times. Therefore, 
possible truncation effects will show up during the time evolution only, but not 
already initially, since the initial ensemble is treated correctly in the Hartree 
approximation. Furthermore, Gaussian ensembles are the ones that are often 
considered in quantum field theory away from equilibrium. It is straightforward 
to construct an initial density matrix that leads to Gaussian correlation functions, 
and vice versa. Finally, from a technical point of view, Gaussian ensembles 
can easily be implemented in both the truncated evolution equations and the 
numerical evolution obtained by sampling initial conditions. 

In the set of Gaussian ensembles, we choose to take the equilibrium distribu- 
tion function of the unperturbed Hamiltonian H Q , where 



H = dx 



= H + V, V = J dx^\ (27) 
Thus, the initial probability distribution is given by 

p[n(x), <j>(x)\ = Z^ 1 exp[-H /T ], Z = J VnV(j) exp[-H /T ], (28) 
where <p{x) = (f>(x, 0), ir(x) = ir(x, 0), and 

j VnV(f) = j YldTr(x)d(f)(x) (29) 

denotes the integral over the initial phase-space. The temperature of the initial 
ensemble is denoted with T . Note that this distribution function is of course not 
the equilibrium distribution for nonzero A. 

Since the ensemble is Gaussian, the only nontrivial correlation functions at 
t — are the two-point functions, and they read 

(0(g)0(g')> = GW(<7, 0)2tr% + q'), G^(q, 0) = T /(q 2 + m 2 ), 

(ir(q)ir(q')) = G n7T (q, 0)2vr5(g + q'), G n „(q, 0) = T . (30) 

Possible variations of this initial ensemble would be to choose different 'initial 
temperatures' T (q) for each momentum mode and the <fi and 7r fields. 

To properly define the model, we formulate it on a lattice in space with spacing 
a. The number of spatial sites is N, such that the volume is L = Na, and we use 
periodic boundary conditions. Due to the finite volume and lattice cutoff, the 
momentum q takes a finite number of discrete values: 

* = {_ + !,. .. >T L (31) 
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and momentum integrals are replaced by sums: / dq/(2n) — > L The dis- 

persion relation is modified due to the Laplacian on the lattice and reads 



-2 , 2 

q + m , 



q 



— (1 — cos aq). 



(32) 



Classical field theory suffers from the Rayleigh- Jeans divergence, which implies 
that the lattice cutoff cannot be taken to zero in a straightforward manner. In- 
deed, the expectation value of the full energy in this ensemble is 



(H) = (H ) + (V)=L 



a 8 



N 



T + a 



3A 



/ 2\2 



(33) 



where the first explicit expression is written such that the extensivity of the 
energy and the linear divergence as a — > are manifest, and the second one on 
the other hand emphasizes 'classical equipartition' in a system with N degrees of 
freedom. In this paper we work at fixed lattice cutoff ma = 0.25 (corresponding 
to a fixed momentum cutoff A = it /a = 4irm). The thermodynamic limit can be 
taken by increasing the number of lattice sites N, keeping the initial temperature 
T fixed. 

The expectation value (<j) 2 ) can be calculated at t — 0: 



1 



^)) = z E g ^°) = tE 



To 



L q q 2 + m 2 



m 



(34) 



The sum is ultraviolet finite. For the lattice sizes we use its value is close to the 
infinite volume value: I = |[1 + a 2 m 2 /4] -1 / 2 ~ 0.5. 

The initial conditions for the evolution equations can now be given explicitly. 
They read 



A(q,0)=G^(q,0) = (q 2 + m 2 )/T , 
B(q,0)=G~l(q,0) = l/T , 



(35) 



and 



C(q,0) = u(q 1 ,q 2 ,q 3 ]0) = v(q 1 ,q 2 ,q 3 ;0) 

= w(q 1 ,q 2 ,q 3 ;0) =y(q 1 ,q 2 ,q 3 ;0) = z(qx,q2,q 3 ;0) = 0. (36) 

The evolution equations are solved using a standard fourth-order Runge-Kutta 
algorithm that is exactly time reversible. 

The full nonlinear evolution is constructed by sampling initial conditions from 
the Gaussian ensemble (j30| ) and solving the equation of motion numerically for 
each initial condition in real space. In order to do this, the action is discretized 
on a lattice in time as well, with step size a < a. The resulting discretized 
equations of motion are of the leap-frog type. We have used a^/a = 0.05 and 
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Figure 2: Time evolution of the two-point function 6^(0, t). We compare the 
Hartree equations (LO, N = 512), the next-to-leading order equations (NLO, 
N = 160), and the results from a Monte Carlo sampling (N = 512, N m = 16,000). 
In Figs. the initial temperature of the Gaussian ensemble is Tq = 5. In all 
figures, the lattice spacing is ma = 0.25. 

0.01, and checked independence of the step size. The ensemble is approximated 
by using many independent initial conditions. The number of ensemble members 
is denoted with N m . For the results presented below we have used N m ~ 2,500 — 
16,000. 

Finally, for the numerical analysis it is convenient to use the mass parameter 
m as the dimensionful scale. Therefore we rescale all dimensionful parameters 
with the appropriate power of m. Rescaled variables will be denoted with a prime. 
Furthermore, the classical equation of motion (^) can be made independent of 
the coupling A by introducing a field cf>' as 



The rescaled (dimensionless) canonical momentum is ir' = (3A/m 4 ) 1/,2 7r and the 
dimensionless energy reads E' = 3XE/m 3 . We define a dimensionless temperature 
T 1 = 3XT/m 3 such that E/T = E'/T'. As a result, besides the lattice parameters 
a' = ma and iV (L' = a'N), only one parameter remains to be specified: 



A larger value of Tq corresponds to a larger effective interaction strength.^ 
4 A quick way to change to primed variables is to put m = 1, A = 1/3. 




(37) 
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Figure 3: Time evolution of (4>' 2 (x,t)) = V 1 Y^qG'^q^i). The initial value is 
(4>' 2 ) = TqJ = 2.5. The time-averaged values in the interval 20 < mt < 50 are 
2.049 (LO), 2.114 (NLO), and 2.126(5) (MC). The Hartree fixed-point value (see 
below) is 2.056. 

5 Early and intermediate times 

We are now fully equipped to compare the time evolution using the different 
methods. For shortness, we refer below to the evolution obtained in the Hartree 
approximation as LO (leading order) and the evolution from the equal-time effec- 
tive action truncated up to four-point couplings as NLO (next-to-leading order). 
Although in principle there is no small parameter governing the truncation, we 
use the labeling common in the large Nj expansion, since the Hartree approx- 
imation is closely related to the leading-order contribution in 1/Nf. The full 
nonlinear evolution of a sample of initial conditions is denoted with MC (Monte 
Carlo). 

The equal-time two-point function at zero momentum G^q = 0,t) is shown 
in Fig. 0, for a typical choice of parameters. The numerical integration of the 
NLO evolution is rather time consuming, due to the presence of the four-point 
couplings u, v, w, y, and z that depend on three independent momentum vari- 
ables. Therefore it is not possible to take as large volumes as in LO and MC. 
However, we have checked that at this stage this does not affect the comparison. 
It is clear that the first few oscillations are well approximated by both LO and 
NLO dynamics. We will refer to this period as the early-time regime. It is visible 
that the Hartree evolution underestimates damping whereas the size of the os- 
cillations in NLO remains comparable with the MC result much longer. Around 
mt = 10 the periods of oscillation in NLO and MC evolution start to disagree. 

The crucial quantity in the mean-field approximation (f|) is the field squared 
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Figure 4: Volume dependence: NLO evolution of (<p' 2 )- In a smaller volume the 
evolution deviates earlier from the large- volume limit. Damping can be seen only 
in the period before the evolution starts to deviate. 

((f) 2 (x,t)) = L^ 1 YtqGffxpfajt) which is presented in Fig. Again we see that the 
first few oscillations are in good agreement. Then a difference becomes visible 
between LO on the one hand, and NLO and MC on the other hand: at LO the 
amplitude of oscillations reduces slower (less damping) and the averaged value 
lies below the other two. A comparison of the time-averaged values between 
20 < mt < 50 shows that the Hartree result differs by a few percent from NLO 
and MC. The time-averaged value in NLO is surprisingly close to the MC result. 

After the initial reduction the size of fluctuations increases again in NLO 
around mt = 30. To investigate this issue, we study the volume dependence of 
the evolution in NLO. The phenomenon of 'resurgence of fluctuations' is presented 
in Fig. H for relatively small volumes. We see that fluctuations become large again 
after the initial reduction and that the evolution becomes undamped. The time 
when this occurs increases with the volume. We stress that this effect is not 
physical, as it is absent in the MC result, nor is it generated by errors in the 
numerical integration. Rather, it is a property of the NLO approximation. Since 
this seems to be an important limitation for the validity of NLO, we analyse the 
thermodynamic limit and the consistency of the NLO evolution on a quantitative 
level. We compare the dynamics for different volume sizes with the largest one 
that is available (N = 160). For an equal-time observable 0(t), we define the 
difference 

AO N (t) = I' ""''-' """ 1 - (39) 

(^160 )av 

The normalization (Oi6o(*))av is the time-averaged value of (Oi 60 (t)) between 
< mt < 20, and is used to set the scale. We denote the time where AO wit) 
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Figure 5: Thermodynamic limit in the NLO evolution: volume dependence of the 
time Tfxt 7v where AOjv(i) exceeds 0.005 (see text), for O = n' 2 (x,t) and 0' 2 (x,t). 
The dashed line is a straight-line fit through the origin. 



exceeds the conservative bound of 0.005 by The dependence of t N on the 
volume is shown in Fig. [5]. It turns out that is sensitive to the details of the 
evolution, due to the oscillating character of AO wit). This sensitivity is indicated 
with error bars. We see that in a larger volume the evolution behaves better for a 
longer time, and that this time increases roughly linearly with the system size. We 
emphasize that the possibility of taking the thermodynamic limit distinguishes 
the comparison performed in this paper with those where quantum mechanics 
was used to test the truncated evolution. For this aspect classical fields are closer 
to quantum fields than quantum mechanics is. 

At a later moment, to which we will refer as t u , fluctuations grow rapidly and 
the numerical evolution becomes uncontrolled. Typically t u is much larger than 
t/v- We have checked that the time where the evolution starts to behave badly 
is not an artefact of the numerical integration. For instance, reducing the step 
size by a factor of 5 does not affect the results. A similar behaviour has been 



noted before in a system of anharmonic oscillators in + 1 dimensions |15|] . The 
breakdown of the evolution equations beyond leading order at large times has 
also been observed for quantum mechanics |TT] . 



Fixed points in the Hartree approximation 

From the viewpoint of thermalization, the most interesting observable is G nn (q, t). 
In an interacting theory the equilibrium value is G^(q) = T for all q. Away from 
equilibrium we define therefore an 'effective temperature' for a momentum mode 
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Figure 6: Time evolution of the effective temperature T'(t) = N l YjqG'^ 1T {q. l t). 
The time-averaged values in the interval 20 < mt < 50 are 5.293 (LO), 5.318 
(NLO), and 5.321(4) (MC). The value at the Hartree fixed point is 5.291. 



T(q,t) = G n7T {q,t), 
and the average temperature over all modes 

T(t) = N- 1 Y,G m r(q,t)=a(<n 2 ). 



(40) 



(41) 



Because of our choice of initial ensemble all momentum modes have initially the 
same (noninteracting) temperature To, so that G n7r (q, 0) is flat in momentum 
space. Due to the nonzero coupling it deviates from being flat immediately after 
t — 0. For later times deviation from 'flatness' of G 7Tn (q,t) is a good measure for 
deviation from thermal equilibrium. 

In Fig. |6] we show the average temperature T' as a function of time. The 
qualitative aspects comparing LO, NLO, and MC are the same as in Fig. [3[ 
Furthermore we see many rapid oscillations with a small amplitude which were 
absent in Fig. |[ The reason is that (0 2 ) is ultraviolet finite and therefore domi- 
nated by the low-momentum modes, while (ir 2 ) is sensitive to all frequencies up 
to the lattice cutoff. 

Fig. ^ gives the impression that the system establishes a new temperature T' ~ 
5.32 ^ Tq = 5 rather quickly. However, in order to be in thermal equilibrium, 
all momentum modes should have the same temperature. In Fig. |7| the effective 
temperature T'(q,t) for three momentum modes is shown. Perhaps surprisingly, 
we see that for each q T'(q, t) oscillates around a different value. The mean values 
appear rather stable and do not seem to approach each other. This resembles the 
nonthermal fixed points discussed in [|TBj for the NLO equations. 
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Figure 7: Time evolution of the mode-temperature T'(q, t) = G' nn (q, t) for modes 
q/m = 0, rr/2, it. The initial value is Tq = 5 for all q. The LO result is not shown 
for clarity. 

In fact, it turns out that this behaviour can be understood with satisfactory 
accuracy already in terms of fixed or stationary points of the evolution equations 
in the Hartree approximation. The relations between the two-point functions 
at a fixed point (denoted with a star) are readily determined from the Hartree 
equations (]13|), and read 

G* m (q)=Q?Cr„{q), G^(q) = 0, (42) 

with 

= < + §*<^>*- ( 43 ) 

The first equation shows an expected relation between tttt and 00 two-point 
functions, the second expression confirms that the time-reflection odd two-point 
function has to vanish at a stationary point. 

The fixed-point structure of the Hartree equations by itself does not yet con- 
strain the allowed fixed-point solutions completely. However, we can supplement 
the set of equations ( f42|) with the nontrivial combinations a 2 {q), given in Eq. 
(|14"D , that are exactly conserved for each momentum mode q independently. At 
a fixed point this gives a third relation 

gimgiM = «~ 2 (<?)- 

We recall that a(q) can be determined from the initial ensemble 
and ([0]) yields the complete fixed-point solution 



NLO 

MC 




(44) 

. Combining ([42j) 

(45) 
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u*a{q)' 



(46) 



Since G*^(q) is identified with the effective temperature for a mode q, the first 
equation shows that at a fixed point the system will generically be nonthermal, 
since the temperature of a mode depends on its momentum, and that the nonther- 
mal 'temperature profile' follows directly from the initial ensemble. The second 
expression leads to a gap equation, after an integration over q, 



2\* 



dq 
2^ 



GlAq) 



2vr u*a(q)' 



(47) 



since the right-hand side depends on (0 2 )* via Lb* . We would like to stress again 
that for an arbitrary initial ensemble the fixed point in the Hartree approximation 
is determined completely and all its properties can be calculated. 

We will now become explicit and specialize to the Gaussian initial ensemble 
we consider in this paper. From the initial expectation values (|30|) one finds 



(4* 



so that the gap equation reads 



i2\* 



Tn 



dq 1 



2tt u*uj q 



(49) 



It is convenient to introduce the dimensionless parameter A = 3A(0 2 )*/(2m 2 ) 
(0 /2 )*/2, so that the above equation can be written as 



V 



. 7T 



A 



27Tv^TA v 2' V 1 + A' 



(50) 



where F(n/2, k) is the (complete) elliptic function of the first kind, and we recall 
that Tq = 3AT /m 3 . This gap equation can be solved numerically. For Tq = 5, 
we find (0 /2 )* = 2A = 2.06, which can be compared with the time-averaged value 
of the LO, NLO, and MC evolution in Fig. |. 

For the effective momentum-dependent temperature at the fixed point we find 



T\q) = GIM = T, 



0J„ 



1 + 



3A(0 2 >* 



1/2 



(51) 



In order to compare this profile with the numerical results, we calculate the time 
average of T'(q,t) between < mi < 50 for the LO, NLO, and MC evolution. 
The result is shown in Fig. |^. We see that the temperature profile emerging 
dynamically in the Hartree approximation is extremely well described by its fixed- 
point shape (pnp. The result from MC turns out to be remarkably close, implying 
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Figure 8: Temperature profile: effective mode-temperature T'(q) versus momen- 
tum q/m after time averaging over the interval < mt < 50, for the LO, NLO, 
and MC evolution. The fourth line is the analytic expression at the fixed point 
of the Hartree equations, with (4>' 2 )* = 2.06 from the fixed-point gap equation. 

that the full nonlinearity only has a small quantitative effect at this stage, and 
that direct scattering is not very important. Also, the NLO profile is close to 
the MC result, showing that the inclusion of momentum-dependent four-point 
functions in the truncation of the effective action improves the agreement with 
the full evolution. Finally, the frequency of oscillation of the individual two-point 
functions G^i (or A, B, and C) is approximately 2cD*, which can be seen in Fig. 
for G nn (q, t). We refer to the stage where the dynamics is well described by the 
Hartree fixed point as the intermediate-time regime. 

At next-to-leading order, the fixed-point structure changes and becomes much 
more complicated. From time-reflection symmetry, it is clear that the two-point 
function C(q, t), and the four-point couplings v(qi, q2,q3) t) and y(qi, q2,Q3] t) have 
to vanish at a fixed point. This implies that also j(q,t) vanishes. However, the 
other four-point couplings cannot be zero at a fixed point, which can be seen 
e.g. from the dynamical equation for v(qi, q2,q3] t) in Appendix fB|. Therefore the 
fixed point is determined by a set of integral equations. Furthermore, the relation 
of the fixed point to the initial ensemble may be rather complicated. 

6 Late times and thermalization 

The results in the previous section show that in the intermediate-time regime 
correlation functions appear quasi- stationary but are not thermal. In particular 
the time-averaged value of G nn (q, t) is well described by the nonthermal profile 
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Figure 9: Time evolution of the zero-mode temperature T'(0, t) (two upper lines) 
from LO (N = 8192) and MC (N = 128, N m = 12,500). The Hartree fixed-point 
value is 28.5. The average temperatures over all modes T'(t) (two lower lines) 
appear as straight lines with the MC result slightly above the LO result. The 
initial temperature of the Gaussian ensemble is Tq = 15. 

flSTD , determined from the fixed point of the Hartree equations. Since the MC 
profile is, up to a small quantitative correction, in agreement with flSTp as well, we 
infer that this (quasi-)fixed point plays a role also in the full nonlinear evolution. 

The fate of the fixed point can be determined by going to longer times. In Fig. 
we show the evolution of the effective temperature of the zero-momentum mode 
for LO and MC. We have chosen a higher value for Tq than before. Also shown 
are the average temperature over all modes, T'(t). As was already mentioned in 
Sec. H, the equilibrium temperature, which we will denote with T', is established 
very early in the evolution, so that the lines presenting T'(t) appear straight. We 
see that at LO the zero mode remains oscillating around approximately 29.2. We 
have calculated the fixed-point value for Tq = 15 and found T'*(q = 0) = 28.5. 
The full nonlinear evolution, on the other hand, shows a decrease towards T'\ the 
approach to thermal equilibrium. We also see that the damping at LO is unrelated 
to the MC result. It is clear that whereas the Hartree approximation describes 
the early and intermediate regimes qualitatively (or even quantitatively), it is 
not able to move away from the fixed point and the approximation breaks down 
completely in the late-time regime. 

Unfortunately, for accessible volume sizes the NLO evolution cannot reach the 
relevant time scales before becoming unreliable (see the discussion around Figs. 
[|, H). At the largest possible times where the evolution could still be followed, 
we have not been able to see a sign of thermalization. 
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Figure 10: Approach to equilibrium: (time-averaged) snapshots of the 
momentum-dependent temperature T'(q, t) for all the modes up to cutoff A/m = 
An. The curves represent a time average over an interval (i, — At, £j), with 
mAt = 1500, and mt x = 1500 (i = 1), 3000 (2), 4500 (3), 6000 (4), 9000 (5), 
and 15000 (6). Line 6 is hardly distinguishable from the dotted horizontal line, 
which is the average temperature from all modes (T" = 24.01). The parameters 
are T' = 20, iV = 256, N m = 2,500 (in Figs. |M| MC only). 



We continue with MC only. As indicated above, a good observable to fol- 
low during the thermalization stage is the temperature profile G n7T (q,t), which 
should become flat (g-independent). The evolution from the fixed-point profile 
at intermediate times to a thermal profile at late times is presented in Fig. ITDL 
We show the time dependence of T'(q, t) for all modes up to the lattice cutoff, 
averaged over an interval mAt = 1500, for six intervals. In the first few intervals, 
the presence of the nonthermal profile is still visible. As time goes on, the profile 
becomes flatter and flatter. In the last interval shown, between mt = 13500 and 
15000, the profile appears g-independent and can hardly be distinguished from 
the straight line, T 1 = 24.01. We see that all momentum modes obtain the same 
temperature T' roughly at the same time. For a detailed investigation on the 
issue of thermalization in this model concerning other correlation functions than 
G nn we refer to our previous paper PU|. The aspects of thermalization we stud- 
ied there are complementary to our findings here (in Ref. |2(J we focused on the 
independence of initial conditions, the long-time behaviour of temporal averages 
in single 'microstates' and other (non- Gaussian) initial ensembles, and the role 
of the thermodynamic limit). 

To determine the time scale for thermalization, we concentrate on the zero 
mode. In Fig. [TT] we show the relaxation of T'(0,t) towards the average temper- 
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Figure 11: Relaxation of the zero-mode temperature T'(0,t) to the average tem- 
perature for three values of the initial temperature Tq. The average temperatures 
T'{t) from all modes appear as straight lines. Also shown are exponential fits, 
explained in the text. The parameters are N = 128, N m = 5,000 each. 



ature T'(t) for three different initial temperatures Tq. The numerical data are 
fitted with an exponential of the form 

T'(0,t) = T' [l + Ke~ t/r ] . (52) 

The fit is performed over the whole time intervalf] The resulting relaxation rate 
1/mr is shown in Fig. |12| versus the equilibrium temperature T', for two system 



sizes, at fixed lattice spacing ma = 0.25. No volume dependence is visible. 

In the remainder of this section, we discuss the thermalization time scale. In 
a quantum theory, the relaxation rate is in general related to the imaginary part 
of the self-energy ||21|| . A recent analysis, applying the relaxation-time approx- 
imation to a weakly coupled quantum scalar field in 3 + 1 dimensions, can be 
found in [23|. At weak coupling, the rate is determined by the imaginary part 



of the setting-sun diagram, describing an on-shell two-to-two scattering process, 
^p + ^k — ► ^q + ^p-k-q- In Appendix |C] we show how to implement a relaxation- 
rate (linear response) approximation for a classical field close to equilibrium. We 
also give the calculation of the imaginary part of the classical setting-sun diagram 
in 1 + 1 dimensions, taken on-shell for arbitrary external spatial momentum. 

One should note, however, that in 1 + 1 dimensions on-shell two-to-two scatter- 
ing is special, since the energy-conservation relation has only two simple solutions: 
q = —k and q = p. Both solutions give u p +ujk — > u} p +Uk, and scattering events of 

5 We also checked for possible power law corrections to the exponential relaxation, but found 
no indication for those. 
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Figure 12: Relaxation rate 1/mr of the zero-mode T'(0,t) versus the final equi- 
librium temperature T' . The initial temperatures are Tq = 2.5, 5, 7.5, . . . , 22.5, 
for A" = 128 using N m = 5,000 initial conditions for each temperature, and 
Tq = 5, 10, 15,20 for iV = 256 (iV m = 2,500). Statistical errors are smaller than 
the symbol sizes. The line is a phenomenological power law fit of the N = 128 
data, 1/mr = 5.8 ■ 10" 6 (T') L39 . 

this type do not change the population numbers for the participating momentum 
modes. We find in Appendix that the naively computed 'relaxation rate' is 
several orders of magnitude bigger than the thermalization rate observed in the 
numerical simulations. 

Processes with nontrivial momentum exchange are necessary for thermaliza- 
tion. This goes beyond two-to-two scattering and occurs only at higher order. 
It is possible that a resummation of the self-energy to two-loop order (i.e. in- 
cluding a resummation of the setting-sun diagram) will determine the relevant 
time scale. This is suggested by the results found in |Tj|, where (quantum) evo- 
lution equations for the two-point function are solved. Those equations contain 
a setting-sun type contribution with fully dressed (in a self-consistent manner) 
propagators. This would also cure the collinear singularity that appears in the 
lowest-order setting-sun diagram, both in the quantum and in the classical theory 
(see Appendix H). 

7 Outlook 

We have investigated the nonequilibrium time evolution of correlation functions 
in field theories. In order to gain an understanding of the validity of approxima- 
tions often used for quantum fields away from equilibrium, we argued that it is 
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useful to consider the equivalent problem in classical field theory. Taking 1+1 
dimensional classical <j) 4 theory, discretized on a lattice, as a simple (but often 
used) model, we have implemented a Hartree approximation and a truncation for 
equal-time correlation functions containing scattering. The truncated dynamics 
was compared with the fully nonlinear results from a numerical sampling of initial 
conditions. We believe that our findings are relevant for the 3 + 1 dimensional 
quantum theory as well. 

For the investigated initial Gaussian ensembles the evolution at early and 
intermediate times is well reproduced by both truncations. We found that the 
Hartree approximation underestimates damping. This difference becomes more 
pronounced at larger coupling. The inclusion of scattering leads to a quantita- 
tively better agreement with the full numerical evolution. 

There is an intermediate-time regime in which correlation functions have non- 
thermal values. The characteristic behaviour in this regime can be understood 
from the presence of a nonthermal quasi-stationary point, which is close to a 
fixed point in the Hartree approximation. As a consequence of (infinitely many) 
conserved correlation functions, the Hartree approximation cannot move away 
from this fixed point and is therefore unable to describe thermalization. 

The quartic approximation to the equal-time effective action (NLO) suffers 
from a different problem. For a finite system (N finite) the initially very successful 
description of the evolution moves away from both the infinite- volume evolution 
at NLO and the numerical results at some time t^- For t > fluctuations grow 
that are related to the truncation and not to the physical system. We find that 
£tv is proportional to N, such that the problem may disappear for iV — > oo. In 
practice this is of little help, since only finite N can be investigated numerically. 
At a time t u 3> the evolution becomes uncontrolled and cannot be followed 
numerically. For accessible values of N, both and t u are found to be (much) 
smaller than the typical relaxation time for thermalization. For this reason, we 
do not know whether NLO is in principle able to describe thermalization or not. 

Since the late-time regime can certainly not be described by the Hartree 
approximation, and the quartic approximation including scattering becomes un- 
reliable on the relevant large time scales, we conclude that the regime of thermal- 
ization is still unsolved in the approach using evolution equations for equal-time 
correlation functions. 

Concerning the late-time regime, promising results for translationally invari- 
ant ensembles have been obtained recently using an unequal-time formulation 
that is time-nonlocal |fL8| . It would be interesting to implement the method em- 
ployed there in a classical theory and carry out a similar comparison as we did in 
this paper for the equal-time formulations. An open question in that respect is 
whether a successful truncation can be found in an equal-time formalism as well 
or if the nonlocality is crucial for an analytical description of the thermalization 
regime. 
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A Symmetry at leading order 



The Hartree equations conserve the combination a(q), defined in Eq. fllil) , for 
each q. This can be understood as follows ||. The effective equations can be 
rewritten by introducing a set of complex variables £(q,t) as 



G^(q,t)=C(q,mq,t), 



1 r: 



In terms of these the evolution equations (|13j) read 



to „ 



sa/^kmi'W*)- 



v q ' 2"J 2tt 

This equation can be derived from an effective Lagrange density 



C " ' 27 



(53) 



(54) 



(55) 



which has a global (i.e. time independent) symmetry £(q,t) — > exp[i6(q)]£(q, t) 
for each g. The corresponding conserved charge reads 



Q(q)=i C(q,mq,t)-C(q,mq,t) 



(56) 



This charge is in fact the conserved quantity and Q(q) = 2a x (g). 



B Evolution of the four-point couplings 

In this Appendix we list the equations that determine the time evolution of the 
four-point couplings. Again, these equations are slightly simpler than the ones 
that can be obtained from setting Nf = 1 in the equations for the 0(Nf) model 

0- 

The evolution equations read: 

d t u(qx, q 2 , q 3 ) = [^(qi, <h, q 3 ) + 4AC(gi) - AXC(q 2 )S 1 (q 1 + q 2 , g 3 )] S YM> 
d t v(qi, 92, qz) = [2a)J 2 tw(gi, q 2 , qz) - 4u(gi, q 2 , qa) + 4AS(gi) 
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-7(?iM?i, <?2, 93) - ^X{B(q 1 )S 1 (q 1 + q 2 , q-s) 
+C(q A )S 3 (q 1 + q 2 , qi) + 2C(q 2 )S 5 (q 2 + q 3 , 9i)}]<jym> 

d t w(q u q 2 , q 3 ) = [3Cj^y(q u q 2 , q 3 ) - v(q u q 2 , q 3 ) - 2v(q 2 , q A , gi) 

-{7(^1) + l(q2)}w(q l7 q 2 , q 3 ) - \{C(q 3 )S 4 (qx, q 2 ) 
+8B(q 2 )Ss(q 2 + q 3 , qi)} - 4A J B(gi)5 3 (-gi - q 3 , qi)] S YMi 

d t y(qi, <?2, ^3) = [4Dj 4 2(gi, q 2 , q 3 ) - 2w(q u q 2 , q 3 ) 

-b(qi) + 7(^2) + 7te)Mgi, 22, ga) - ^(q^s^, q 2 )\ SYM 

d t z(q u ?2, qs) = ~[y(qi, <?2, 9a) + 47(gi)^(gi, 92, Qs)]sym^ 

where the subscript sym implies symmetrization with respect to the appropriate 
permutations of q%, q 2 , q 3 , and g 4 = —{q\ + q 2 + q 3 ). Here we have used the 
one-loop integrals 



3 r 

Si(qi,q 2 ) = - 2TT5(p 2 -pi - qi) G H (pi)G H {p 2 ) [6u(p 2 , -pi, q 2 ) 

4 Jpi,P2 



-3c(px)i;(-pi,p2,?2) + c(pi)c{p 2 )w(-p 1 ,p 2 ,q 2 )}, 
S 3 {qi,q 2 ) = -; / 2tc5(p 2 - p 1 + q{) G H (p l )G H (p 2 )[lv(q 2 , -pi,p 2 ) 

4 Jpi,P2 

-8c(p 2 )w(p 2 , q 2 , -pi) + 7c(pi)c(p 2 )y(p 2 , -p u q 2 )}, 
S4{qi,q 2 ) = 3 / 2tc5{p 2 -pi + q 1 + q 2 ) G <t>tt> (pi)G H {p 2 )[w{qi,q 2 ,p 2 ) 

J Pi ,P2 

-3c(pi)y(gi, g 2 , -pi) + 6c(pi)c(p 2 )2(P2, -Pi, <?i)], 
Ss(qi,q 2 ) = - 2ir5(p 2 +p 1 + q 1 )G (f , (l) (p 1 )G 4>lj ,(p 2 )[v(q 2 ,-pi,-p 2 ) 

4 Jp\,P2 

-2c(p t )w(-pi, 92, -pa) + c(pi)c(-p 2 )|/(g2, -Pi, -Jfc)]. 

C Relaxation-time approximation 

For a quantum field slightly away from equilibrium the relation between the imag- 
inary part of the self-energy and the thermalization rate for the single-particle 



distribution function has been pointed out long ago by Weldon [^TJ. A more 



recent detailed analysis can be found in [^] for a scalar field in 3 + 1 dimensions. 



Here we briefly outline the arguments of |22]| in d + 1 dimensions and then show 
how to adapt those for the classical theory we consider here. This analysis is 
valid for a weakly coupled plasma, close to equilibrium. 

A dynamical equation for the (quasi-)particle distribution n(jp,t), describing 
the relaxation towards the Bose distribution 7Zb(<^ p ) = l/[exp(u p /T) — 1], can 
be obtained perturbatively, using the Heisenberg equations of motion and resum- 
ming hard thermal loops if necessary. In the relaxation-time approximation, the 
distribution function for a momentum mode p is written as 

n(p,t) = n B (u p ) + 8n(p,t). (57) 
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All other modes q ^ p are assumed to be in equilibrium. Skipping many steps 

(58) 

ImS(w p ,p) 



22j] , the resulting linearized equation is 

d t 5n(p,t) = -T(p)6n(p,t), 



with the relaxation rate 



t(p) 



(59) 



This rate is twice the plasmon damping rate 7(p). For a weakly coupled scalar 
field with a 4 interaction in 3+1 dimensions, T(p) is determined by the imaginary 
part of the setting-sun diagram. 

For the classical theory our interest is in the evolution of the momentum- 
dependent 'temperature' T(p, t) = G^p,^. We may obtain this correlation 
function from the classical unequal-time two-point function S^x — y;ti, £2) — 
(0(x, ti)<f)(y, t 2 )) by a spatial Fourier transform as 



T(p,t) =d tl d t2 S{p;t 1 ,t 2 



t!=t 2 =t 



(60) 



A calculation of S(p; t±, t 2 ) using classical perturbation theory to second order in 
the coupling constant can be found in [^] for d = 3 (see also [24|). The resulting 
expressions are similar to the quantum ones. The most important change is that 
the Bose distributions are replaced by classical distribution functions 



T 



0J r 



(61) 



The relevant second-order contribution reads [p3 
Sa(p;*i,*2) = 



poo 

- / dt'dt" G R (p, h - t')Y, RA {p, t' - t")S (p, t" - t 2 ) 
Jo 

POO 

- / dt'dt" S (p, h - t')E Acl (p, t! - t")Gtf (p, t" - t 2 ) 
Jo 

poo 

- / dt'dt" G H ( P , ti - t')E FA (p, t! - t")G$(p, t" - t 2 ) 
Jo 



where 



S (p,t) = n cl (uj^ 



is the free two-point function. The free retarded Green function reads 

G^p,t) = e(tf^^ = G^p,-t), 

Up 

and the self-energy corrections are, in d + 1 dimensions, 
£h, c i(p,£) 
£i?, c i(p,t) 



(62) 



(63) 



9A 2 


/ 
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A,q 


9A 2 


/ 


6 


Vk,q 
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with 

d d k . , 

Wf- < 64 > 

After performing all the time integrals and taking the derivatives as in Eq. (|S0D , 
one may follow the arguments given in [22] for the quantum theory to find the 
evolution for ST(jp, t) = T(p, t) — T in the relaxation-time approximation. The 
result is 

d t 8T{p,t) = -T cl {p)5T{p,t), (65) 

with 

r = _ Im^K, P ) (66) 
On-shell, only the two-to-two scattering contribution in the retarded self-energy 



is kinematically allowed. This can be written as |25|, |26 



ImS RiC i(u;p,p) = 

J rf$i 23 (p) 2n5(u p + u k - cu q - ^ r )n c i(^ k )n cl (^ q )n c i(u; r ), 



9A 2 u p 



where 



^(p) = Jl ( otl " k " q " ^ < 67 ) 

(27r) a 2a;k (27r) a 2cG> q (27r) <:1 2ci; r 

Specializing to d — 1, we find 

9A 2 T 2 

Tci(p) = -^HP), (68) 

with 

/(p) = r dk r <i q s ^ +u *- u r u '- k - } - m 

J-co J~co U k U q U p _ k _ q 

Note that the momentum integrals are ultraviolet finite. The integral is invariant 
under p — > — p, so we may restrict ourselves to p > 0. The g-integral can be 
performed using the delta function, which has support at two separated points 
only, q = p and q = —k. The result is 

2 roo 1 1 
J(p) = — dk- 



oo u k | A;c<jp + puk | 

U p J-oo U)% W%-U* U p J-p U% Ul-U* 

The remaining integrals contain a collinear singularity when k — ► — p, leading to 
a logarithmic divergence. This singularity is present in the quantum self-energy 
as well. We regulate this in an ad-hoc manner by modifying the integration 
boundaries to — p ± /i, with \l <C m a small cutoff. 
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The integrals are straightforward using partial fractioning, and the final result 

is 



9A 2 T 2 1 

r c i(p) 



lQnm 2 Up 



P P , 
— 1 H arctan h In 



(71) 



The rate is positive for sufficiently small fi (fi/m < 0.7). Two limiting cases are 

„ , , 9A 2 T 2 2m „ , , 9A 2 T 2 

r c i = - In — , r d (p » m) = 3 2 . (72) 

Up to now we have silently neglected the corrections to the mass parameter due 
to interactions. If we denote the resummed mass parameter with M, m has to 
replaced by M in the expressions above. 

Let us conclude by noting that loosely identifying fi with T c i itself gives a 
relaxation rate r c i(0) which is at least two orders of magnitude too big, when 
compared to the numerical results for the thermalization rate 1/r. 
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